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Abstract 



^^ ■ Once superpartners are discovered at colliders, the next challenge will be to 

C^^ ■ determine the parameters of the supersymmetric Lagrangian. We illustrate 

Q ' how the relative phases of the gluino, SU{2), and U{1) gauginos and the 



Higgsino mass parameter fj, can be measured at a hadron collider without ad 



^i I hoc assumptions about the underlying physics, focusing on Fermilab. We also 



discuss how the gluino and LSP masses can be measured. 
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I. INTRODUCTION 

Presently there is considerable indirect evidence that the Standard Model (SM) is ex- 
tended to a supersymmetric SM (SSM), and that at least some superpartners are light 
enough to be produced at the LEP or Fermilab colliders. There is some motivation that 
gluinos are light and will be produced copiously at Fermilab in the next run [jl|, based on 
the success of supersymmetry in explaining electroweak symmetry breaking, which relates 
superpartner masses to the Z boson mass. Once superpartners are discovered, the next step 
will be to measure the parameters of the Lagrangian, and then to learn the underlying theory 
that leads to such a Lagrangian (technically, the softly broken supersymmetric Lagrangian). 

Before using existing information to constrain the Lagrangian, the minimal case has 
over 100 free parameters [0, of which at most two are directly measurable, the gluino and 
gravitino masses, since these are the only ones that do not mix to form mass eigenstates. 
The parameters consist of (at least) 33 masses, 40 phases, and 32 super-CKM type angles 
relating flavor and mass eigenstates. The phases can all induce CP violating effects, and it 
has recently been argued that all CP violation could be explained by these SUSY soft phases 
. The phases affect much more than CP violating observables - cross sections, branching 
fractions, radiative corrections, the LSP relic density, LEP Higgs boson hmits, etc. Q-Q. 
In general, superpartner production cross sections, decay rates and kinematic distributions 
depend on the phases, as we illustrate in this paper. The same gluino phase 03 that we 
consider is also measurable in CP violation in the Kaon system [^]. 

In the past |0], it was argued that the supersymmetric soft phases must be small because 
they would otherwise induce neutron and electron electric dipole moments larger than the 
experimental upper limits . However, recently it has been shown |]8HT0| that large phases 



are not disallowed by current data provided that reasonable relations exist among them. 
In fact, one expects relations among parameters in typical theories. For example, a model 
produced by embedding the MSSM on a particular configuration of D-branes naturally has 



large phases and small EDMs [|TT|. The phases of the supersymmetry Lagrangian may or 
may not be small - it will be necessary to measure them to find out. 

In this paper we illustrate how the presence of phases affect observables and can be 
measured in gluino production and decay. The gluino production cross section does not 
depend on the phases, but the gluino decay does. To understand the essence of what happens, 
we consider the tree level production of gluino pairs followed by the decay g -^ ggNi, 
where Ni is the lightest neutralino and the LSP. The generalization to other production and 
decay channels is straightforward but complicated. For simplicity, we assume that squarks 
are significantly heavier than gluinos, so that squarks decouple from the production cross 
section. 

The phases of interest enter initially in the gaugino and higgsino mass parameters M^e^'^^, 



M2e*'^^, Ml 6**^% and /xe*''^^ in the Lagrangian. For our example, the gluino phase enters 
directly, and the others through the neutralino mass matrix 



Y 



/ Mie^-^i -Mzs/ssw MzCpSw \ 

M2e*^2 Mzcpcw -MzCpCw 

-Mzspsw Mzspcw -/ie*"^^ 

V Mzcpsw -Mzcpcw -/xe^'^" ) 



:i-i) 



where our notation is sp = sin/3,c^ = cos/3, svk = sin^vK and cw = cos6'vi^. This matrix 
can be diagonalized by a 4x4 unitary matrix N, i.e. Yd = N''"YN. U(l) symmetries of 
the Lagrangian allow reparameterizations 0,^,13], so that only the combinations 03 — 02, 
03 ~ 01; and 03 + 0^ are invariant and, thus, observable. 

For the case considered here, we will see that only one phase can be measured, 0eff, 
which depends on all of these Lagrangian phases (and also on tan/3 and masses). When 
additional decay channels are included, different combinations of the Lagrangian phases can 
be measured. Measurements at Fermilab alone could establish that (pes was non-zero, and 
therefore that at least one phase was non-zero. That would mean that at least one phase 
was fairly large in order to give a significant numerical effect. It is more difficult to show 
that the phases are small, since some combinations of them may be small but not others. 
A combination of measurements from electric dipole moments, 6-factories, and high energy 
colliders producing superpartners would establish that the phases were small if they indeed 
were. 

It is important to emphasize that we are following a general procedure, with no model 
dependent assumptions about phases being small or related to one another, or masses being 
degenerate, etc. Our simplifying assumptions are used only to present the results simply; the 
procedure can be followed without simplifying assumptions once real data once is available. 

In addition to phases, two quantities are particularly important to measure because they 
affect many other results. One is the LSP mass. The LSP is a very good candidate for the 
cold dark matter of the universe. If superpartner events are observed with an escaping LSP 
it will, of course, greatly encourage us to indeed believe the LSP is the CDM. But that is 
not established until a calculation of the relic density is carried out that shows ^lsp ~ 1/4, 



and the calculation cannot be done |T^ until the LSP mass and couplings are measured 
and perhaps some superpartner masses. If the LSP candidate interactions are observed 
in explicit dark matter detection experiments, the LSP mass deduced from those can be 
compared to the LSP collider mass. In order to measure the LSP mass we provide a two 
step procedure. 

The second important quantity is tan/?, but it is very difficult to measure. Here, it 
occurs in the A^i couplings (see next section) to quark-antiquark pairs, which enter into 0eff 
along with phases and masses. To the best of our knowledge, tan/3 can only be definitively 



measured at a lepton collider with a polarized beam and sufficient energy to produce some 
superpartners 0]. 

In the next section we explain the formalism and derive the relevant distributions. 

II. RELEVANT LAGRANGIAN AND ANALYTICAL RESULTS 

The general gluino mass term in the MSSM Lagrangian is 

C D -i(M3e*'^^A,A3 + Mae-^^^A^A,), (2.1) 

where A^ is a Majorana spinor. With the field redefinition il}g = GXg, G = e*~ , the kinematic 
mass of the gluino is real. After this redefinition, the quark-squark-gluino vertex is 

C,gg = -g,V2T^,{G*^^PLqMl ' G^Pnqfqtn + (MPLQagl - G4Pb3^, (2.2) 

where the lower index i is the fiavor label, and the upper indices j and k label the color of 
(s) quarks. 

The quark-antiquark-neutralino vertex is 

^.HN. = -gV2fl^liPRNiiL + 9V2fP\PLN,qiR + h.c, (2.3) 

where 

//- = n^Nu - tan^^(T3, - e,)Nn /f = tan^vye^iVn. (2.4) 

The Nij are the complex elements of the unitary matrix which diagonalizes the neutralino 
mass matrix, and depend on the phases 4>i,4>2, (p^i- 

The production of gluino pairs does not depend upon the phases, and can be calculated 



using the standard CP-conserving Feynman rules of the MSSM |Tj]. In the limit that the 
squark masses are very heavy, the spin structure is analogous to that of top quark pair 
production [|16|. The details of this calculation are contained in the appendix. We have 



included this spin structure in our calculation to test whether it infiuences any physical 
observable. In the end, we observed no spin correlations, so the calculation is only mildly 
interesting. 

Apart from an overall color factor, the polarized decay amplitude squared for a gluino 
with spin vector s is: 

\9{s) ^qiQjNil'^ = \0\l + 26L\mg{\aL\'^pi ■ PjS ■ Pi - |/3l|V ■ PiS ■ Pj) 

+2miRe{aLpl)iPg ■ PjS ■ Pi - Pg ■ PiS ■ pj) 
-2m^lui{aLl3l)e^'""'Pg^,s,p,pP^^^ + (L -^ R). (2.5) 



The spin independent piece |Op = |0|| + |0||.is 

- (^laL^i- m\){m\ - t) + \^l\^{u - vr?^{m\ - u) 

+2Re(aL/32)mim^s') + (L -^ i?), (2.6) 

with the couphngs 

ccL = 2g,gGfl^*/Ml an = 2g,gG*f^/Ml 

5l = 1,5r = -1 PL = al l3R = a*R, (2.7) 

and the kinematic variables s = (Pi + Pj)^, i = {Pg—Pj^, and u = {pg—PiY- The gluino and 
neutrahno four- momenta are denoted by Pg and pi, respectively, with p^ = m? and p\ = m\. 
Pi and Pj denote the quark and antiquark four- momenta, and we assume pf = p] = Q- These 
expressions have been summed over the spins of the final state particles. In the limit that 
the phases are zero, these expressions reproduce the standard results. A similar formulae 
was derived earlier for the three-body decay of a heavy neutralino into a lepton-antilepton 
pair and a lighter neutralino 0. In writing out the expressions (|2.7|), we have assumed 
universality of the squark masses (except for the top squarks, which are not important for 
the analysis considered here) 



"n^iiL = ML,mg^^ = M, 



R 



where index i runs through five flavors of squarks, and that the squark masses are much 
larger than any momentum g^ exchanged in the decay, making the replacement 



2\— 1 , s^— 2 



(g^ - M^' -^ -M 

Our methodology is valid regardless of this assumption, and it is not necessary for analyzing 
the real data. To generalize these expressions, one would simply replace M| ^ by M| ^ — i 
in aL,R or by M|^^ - m in Pl,r- 

Eq. ( p.6| ) can be simplified, using the fact that |q;l| = \Pl\ and \aji\ = \Pr\ (which is only 
true with our simphfications) , to 



-(IttLp -F Ittijp) {t-ml){ml-t) + {u- ml) {ml- u) +4 COS (peamiTJigS . (2.8) 

The effective phase (p^s is defined by 



\aL\^ coscpL + laRl cos (pR 
cos0eflf = , — f^— n 13 (2-9) 

where (pL is arg(Q;2,) and (J)r is arg(a;ij). There is a summation over the index i, which 
specifies the quantum numbers of each squark flavor. 



The expressions for the production of gluino pairs and the decay outhned above are the 
starting point for our phenomenological analysis. Focusing on the decays, we note that 
the effect of phases appears in Eq. ( [^.8D as the coefficient of the quark-antiquark invariant 
mass. Therefore, we expect to observe sensitivity to phases in that distribution. A similar 
observation was made for the case of three-body decays of neutralinos [Q. The form of 
Eq. ( ^.9[ ) guarantees that phase effects will not decouple if one of the squarks is much 
heavier than the other. We need not have Af^ ~ M^ to obtain a sizable effect. 

Although most generally there should be four relevant phases, by suitable field redefini- 
tion, it can be shown that only combinations 0i — 03, 02 " 03 and 03 + 0^ are relevant for 



this process. In fact, these are all R-invariant physical observables P,|9|JT3[], and therefore 
cannot be further rotated away by field redefinition. 

Finally, the phases are a manifestation of CP-violation, and their presence should be 
manifest in a CP-violating observable. One example considered here is the contraction of 
particle four-momenta with the 4-dimensional Levi-Cevita tensor. 

III. PHENOMENOLOGY AND MEASURING PARAMETERS 

We concentrate on the signature of 4 jet plus missing transverse energy If^x which is 
naturally expected for the process considered here. To define the signal, we specify the 
following cuts, motivated by the D0 Run I multijet and I^t analysis [|TH|: 



^T > 75 GeV.KT > 100 GeV, |r/J"*| < 2.5, p^^* > 15 GeV, Ai^^^ > 0.5, 

0-1 < ^<^je,^. < - - 0-1' J(^ - ^K4r^' + ^'^Ij, > 0-5 (3-1) 



where ARij is the rj — (f) separation of jets i and j, Tir is the scalar sum of j?^ excluding the 
leading jet in Et, and 1 and 2 are subscripts denoting the leading and next-to-leading jet 
in Et- The A0 cuts are necessary to reduce the fake IpT backgrounds that arise when part 
of a jet is lost. For this analysis, we ignore the fact that some background events will pass 
these cuts. Their effect can only be subtracted statistically, and they will degrade the results 
presented here. The application of cuts biases the final event sample, and the details of the 
proposed measurements presented here will depend upon them. The exact cuts needed to 
reduce backgrounds (primarily the cut on If^x) will not be known until the experiments begin 
to collect data at high luminosity. These cuts are meant to be representative of what will 
actually be done. Finally, we smear jet energies by a Gaussian resolution [aE = SOy/E, E 
measured in GeV) typical of Tevatron experiments. 

Some of the measurements will benefit if the flavor of the jets can be tagged. If all of 
the jets originate from light quarks, then there is a threefold ambiguity in assigning pairs 
of daughter jets to a mother gluino. We have tested various kinematic quantities, and find 



that choosing the combination which minimizes the sum of the square of the jet pair masses 
does fairly well in picking the correct one. However, there is a substantial degradation in 
the measurements with this choice. On the other hand, if two of the jets are fe-tagged, then 
the ambiguity is substantially reduced. We estimate that this can retain roughly a fraction 
of 2 X .2 X .8 X .36 = .12 of all events, where .36 is an estimate of the double 6-tag efficiency. 
Charm quark tagging is also possible. Of course, heavy flavor decays have to be corrected 
for the energy lost to neutrinos, an effect which we do not simulate here. Finally, jet charge 
can be measured for even light quarks, using the fact that a u-quark will fragment to a 
leading 7r°, vr"*", i^"*" with (very) roughly the same probability, while a (i-quark will fragment 
to a leading vr"*", vr", i^°. Therefore, if four jets are observed with leading tracks of charge 
+, — , and 0, the combinations (+, — ) and (0, 0) would be preferred. 

As explained earlier, we have made the simplification that the squark masses are universal 
and heavy in order to present the results simply. Since the mass splitting AM-^^ = rrig—mff^ 
can be measured in a straightforward manner, we choose the value 145 GeV for our numerical 
work, and consider gluino masses from 200 to 350 GeV in steps of 25 GeV. The final results 
will not depend dramatically on the value of the mass splitting, provided that it is large 
enough so that jets from the gluino decay are measurable. 

A. Coupling Measurements 

Even with the simplifications made in this analysis, the gluino decay distribution depends 
on several parameters of the MSSM Lagrangian. From Eq. ( p.8|) , it is clear that we are 
sensitive to only an overall normalization and the relative strength of the s term to i and u. 
We can perhaps measure the overall normalization through the effect of the decay width on 
the shape of the gluino invariant mass distribution, but the width is much smaller than the 
typical energy resolution. The branching ratio for the qqNi decay may be inferred from the 
number of events observed, but one would have to measure other decay modes to make use 
of this information. Therefore, we concentrate on methods to determine 0eff- 

One method is illustrated in Figs. |T] and ^ In Fig. |T], we plot the invariant mass dis- 
tribution of the jets assuming that the jet pairs can be unambiguously identified with a 
particular gluino. This would require, for example, that two jets have a heavy fiavor tag and 
two jets have an anti-tag. The three curves (solid, long-dash, short-dash) show the effect 
of (pes = 0, 7r/2, 77. From 0eflF = ^ vr, the peak has shifted by approximately 25 GeV. If 
limited statistics require a combination of all data, so that tagging cannot be used, the one 
method for pairing the jets is to minimize the sum of the squares of the invariant masses. 
The resultant invariant mass distribution is shown in Fig. Q. There is still a phase depen- 
dence, but the miscombination has washed the effect out somewhat (the peak has shifted 
roughly 10 GeV). 



To estimate how well we can actually distinguish these distributions, we generated sets of 
fake data based on the distribution for 0cff = tt and compare them to the ideal distribution 
for 0eflf = 0. The distribution of x^ based on this comparison gives some indication of the 
sensitivity of the experiments. We quote the value of x^ that contains 5% of the total area 
of the x^ distribution. We could also quote the median value of the distribution, but this 
would give no indication as to the size of the low-side tail. The interpretation of the x^ 
depends upon how the fake data was binned in a histogram. For this case, x^ = 1.9 and 
2.5 corresponds to a 95% and a 99% confidence level, respectively. Therefore, a x^ that 
is far from 2.5 indicates a good separation of the two distributions. A x"^ near 1 means 
the two distributions are consistent with each other. We first consider the case when the 
two daughter jets can be paired correctly with the mother gluino. Assuming an efficiency 
of 0.2 for correctly tagging the jet pair, nig = 250(350) GeV, and 2 fb~^ of data, the 
comparison yields x^ = 11.7(1.3). The value x^ = H-^ for rrig = 250 GeV shows that the 
two distributions are quite distinct with just 2 fb^^ of data. On the other hand, with only 
2 fb~^ of data, the distribution for ^efr = vr and rrig = 350 GeV is fairly consistent with 
that for 0efT = 0. However, with 10 fb"^ of data, x^ = 6.7. Accepting miscombinations of 
the jet pairs, but an efficiency of 1.0, the 2 fb~^ numbers are 23.7(2.1) for rrig = 250(350) 
GeV. With 10 fb~^, the x"^ for f^g = 350 GeV increases to 11.0. In all cases, we assume a 
branching ratio of 1.0 for the four jet decay mode. In general, we find that the increase in 
sensitivity scales linearly with integrated luminosity. Despite the fact that miscombination 
of the jet pairs distorts the invariant mass distribution, a large (pes dependence still remains, 
so that the additional statistics makes for a better separation of the two hypotheses (0eff = 
or vr). 

From pure kinematic considerations, the invariant mass distribution of the jet pairs has 
an end point at AM-^^ = rrig — mi (the end point is smeared if we choose the wrong 
pairing of jets and by energy resolution). For the measurements discussed above, and others 
discussed below, it is important to measure the endpoint with sufficient accuracy. This 
requires a detailed knowledge of both the number of background events and the shape of the 
background distribution (as well as the same quantities for the signal). The cuts of Eq. ( p.l|) 
are somewhat more restrictive than those used in the D0Run I analysis [0, but we can use 
them to estimate the background at 0.5 pb. This means that the signal to background ratio 
for rrig = 250 GeV is roughly 1, decreasing to 1/10 for rrig = 350 GeV. Tighter cuts may 
be desired to establish a signal for heavier gluino masses and to measure parameters. For 
two very different choices of the background shape, and for cos 0cff — 1 and rrig = 250 GeV, 
the endpoint of the invariant mass distribution can be measured to a few GeV. There is a 
systematic shift in the fit endpoint and roughly twice error for when cos0eff — ~1; but this 
is an artifact of approximating the mass spectrum by a straight line near the endpoint. This 
will not occur in a real analysis that includes the full shape of the signal distribution. While 
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the full invariant mass distribution, including miscombinations, gave a better separation of 
different parameters, a tagged sample of events will probably yield a cleaner measurement 
of the endpoint of the invariant mass distribution. 

One distribution that does not depend on the jet pairing or the measurement of the 
endpoint is based on contracting the jet 4-momenta with the Levi-Cevita tensor. We define: 

^'^ u P a 
_ ^pi^p(7PlP2P3P4 fn 2) 

E1E2E3E4 

where the momentum components are measured in the laboratory frame. This distribution 
is shown in Fig. ^. The distinguishing feature is the half width at half maximum for the 
curves, which varies from 1.0 to 1.5 as 0eflf = —>■ tt (this is for the case of rrig = 250 GeV, 
but the variation for other masses is small). We estimate the sensitivity as for the invariant 
mass distribution. For rrig = 250 GeV, x^ = 9.3 with only 2 fb~^, while for rrig = 350 GeV, 
X^ = 1.0(4.8) for 2 (10) fb~^. Therefore, the untagged jet-jet invariant mass distribution 
yields a superior separation of (pcS = or vr compared to any other method considered here. 
The best result will come from combining several observables. 

For Figs. |l]-|], we have used the particular example of rrig = 250 GeV. The invariant mass 
and e distributions discussed above are not very sensitive to this number, so the existence 
of CP-violating phases can be established without knowing the masses in the problem (the 
mass splitting is inferred from the endpoint of invariant mass distribution). In the next 
subsection, however, we will demonstrate that the masses can also be known. 

B. Mass Measurements 

From the observed mass splitting AM-^^ , a measurement of rrig will accordingly give us 
a measurement of mi. Even in the case when the mother gluino is not tagged, comparison 
of different rriij distributions can yield a measurement of AM-^^. This is demonstrated in 
Figs. and ^ previously mentioned. 

By comparing the total number of events to the theoretical prediction for the cross sec- 
tion, one can estimate rUg, and therefore determine all the masses in the problem. However, 
the total number of events really only measures a x BR^{g -^ qqNi), and some decay modes 
may not be measurable at a hadron collider. What is preferable is a statistical measure of 
the gluino mass based on kinematics. As a motivation, we direct the reader to the measure- 
ments of the top quark mass in dilepton events, which can be performed even though there 
are two neutrinos per event [0. Given all the kinematic constraints for ti production and 
decay, there is actually only one unknown, and the momenta of the neutrinos can be recon- 
structed (there are four solutions) for an assumed top quark mass. There is a wide range of 
masses rrit that yield "reasonable" solutions, but a probability can be assigned to each rrit 
by comparing the resultant kinematic distributions to those expected for the assumed rut- 



The case of gluino decay is more complicated, since there are no Mw constraints. Even 
after specifying mg, the longitudinal momentum pz of the LSP's is not fixed. We proceed 
by making simple fits (using our parton level generator) to Pz as a function of nig (the sum 
of a wide and a narrow Gaussian), properly normalized to unity. The individual Pz values 
are not noticeably correlated, and we assume they are independent in making the fit. For 
each event, a range of values for mg is assumed, and an integration is performed over the 
two Pz distributions (the integral is replaced by a finite sum). For each fixed value oip\, pi, 
and rrig, the full kinematics (up to multiple solutions) are reconstructed. Finally, a Gaussian 
weight or probability is assigned to each rrig based on the difference {y/I — 2mg)/{20 GeV). 

The resultant distributions as a function of rrig are displayed in Fig. ^ for the correct 
assignment of jets to a mother gluino (this assumes picking those events which satisfy a 
tagging as outlined above). The same distributions using the minimum m^ algorithm are 
shown in Fig. ^. First, we note that we always reconstruct the correct mass (the peak of the 
weight function) within about 10 GeV for the correct assignment of jets. The case of (pes = tt 
seems to yield a larger mismeasurement, but this is not quite correct. For simplicity, we have 
used the expected pz distribution for the case of (peS = in all cases. However, we showed 
above that the existence of a phase can be inferred independently of the gluino mass, so 
the correct distribution can and should be used for each phase. Secondly, the reconstructed 
value considering combinatoric ambiguities (Fig. ^ is always roughly 25 GeV higher than 
the real mass, which can be corrected. Also, the distributions are much wider, so there will 
be a greater uncertainty in the measured mass. To estimate how well the masses can be 
measured, we fit a Gaussian to the distribution for rrig = 250 GeV around the peak and 
recorded the error on the fit mean. Allowing for miscombinations of the jet pair, the fit 
mean was 272 ± 5(2) GeV for 2 (10) fb~^ of data. Assuming the correct jet pairing, but an 
efficiency of 0.2, the fit mean was 259 ± 9(4) GeV for 2 (10) fb"^ of data. 

The method outlined here is merely a simple demonstration of principle. While details 
such as soft gluon radiation, non-Gaussian errors in energy measurements, and background 
contamination would degrade our results, it is also possible that more sophisticated pattern 
recognition techniques would yield a net improvement. 

IV. CONCLUSION 

In this paper, we have considered some of the effects of CP-violating phases from the 
supersymmetry Lagrangian on the phenomenology of gluino pair production and decay. Such 
phases also affect non-CP-violating observables. Our specific results are for the Fermilab 
collider, but the general approach apphes equally well to the future LHC collider. We 
have also considered the issue of measuring sparticle masses, which is not particular to 
the existence of the phases. Although our analysis uses certain simplifications in order to 
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demonstrate analytically how the phases enter, the method we present here shows that those 
measurements can be done if we have experimental data for that process. Our method can 
be used as a prototype for more complete analyses to treat real experimental data. 

We concentrated on the case of gluino pair production, followed by the gluino decay to a 
quark-antiquark pair and a neutralino LSP. We found that there are simple measurements 
of the gluino decay products which are sensitive to a combination of phases and couplings 
which can be parametrized by an angle 0eff- If 4'es is non-zero, then it is a sign of CP 
violation originating in the supersymmetry soft-breaking Lagrangian. 

We have also demonstrated a method to measure the LSP mass. This is based on first 
measuring the gluino mass rrig, even though there are two invisible LSP's in each event. 
The method is similar to that used to measure the top quark mass in dilepton events: by 
comparing the observed kinematics of an event to the theoretical prediction for a range of 
assumed gluino masses, one can determine which gluino mass is most consistent with the 
data. In our analysis, the most consistent gluino mass is very close to the real gluino mass 
used to generate events. Then, once rrig is known, the LSP mass is determined from the 
endpoint of the jet-jet invariant mass. The measured gluino and LSP mass are not the exact 
quantities that appear in the MSSM Lagrangian. The physical gluino mass rrig can be related 
to the gluino mass parameter M3 through QCD, but this does require some knowledge of 
the squark masses. More measurements will be needed to relate mi to the parameters Mi, 
M2, fi, and tan/3. 

If Supersymmetry is responsible for stabilizing the hierarchy between the weak and 
Planck scales, then there is a strong possibility that sparticles will be produced in the 
near future at the Tevatron. If the gluino is light enough, it will be produced in quantity, 
allowing its discovery and, as argued in this paper, a measurement of its and other sparticle 
properties. Based on our present knowledge, these measurements might indicate that Su- 
persymmetry is responsible for all of CP-violation and/or that a neutralino LSP constitutes 
the dark matter of the universe. 
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APPENDIX: CALCULATION OF PRODUCTION RATES 

To maximize the sensitivity to couplings, we calculated the tree level gluino pair pro- 
duction processes including the spin correlations. We followed closely the calculations of ti 



production in the "off-diagonal basis" [p!5HT7||. When the squark masses are heavy, one finds 
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the same topology of Feynman diagrams for the qq -^ gg subprocess as for ti production, 
and the latter results can be used directly with a simple modification of the color factor. 
The relevant expressions can be found in the above references. To summarize, in the off- 
diagonal basis, when the gluinos are produced in the narrow width approximation, only the 
spin combinations (T, i) or (i, t) occur, where the gluino spin is quantized along the special 
axis that defines the basis (in the rest frame of the gluino, the spin vector makes an angle 
TT — ^ with respect to the gluino direction of motion in the gluino pair rest frame; The angle 
^ is related to the scattering angle 6* in the gluino pair rest frame and the gluino velocity 
via the relation tan^ = y/1 — P'^tanO*). Once the gluino decay is included, there can be 
interference between the (T, i) and (i, T) amplitudes, but these are suppressed by a factor 
of /?^. To simplify our phenomenological analysis, we ignored this interference effect. This 
approximation was justified in the end, since we observed no spin correlations between the 
two gluino decays. Each gluino was decayed with a distribution given by the formulae of 
Sec. II. 

The subprocess gg — > gg has a more complicated spin structure, and is numerically 
less important than the qq subprocess at the Tevatron for the masses considered here. For 
completeness, we included this subprocess in our final results. The structure is almost 
identical to that for gg — > ti, except for the overall color factor and the replacement (9 + 
7 cos^ 9*) -^ (3 + cos^ 9*) in the formulae of Ref. |T^. We included the different contributions 



from (I, t), (I, I), (t, I) and (|, |) and added the decays ignoring the possible interference. 
Again, this is justified by the fact that none of our observables depended significantly on 
the spin correlations. 

For the full calculation, we used CTEQ3M structure functions, NLO running for a^, and 
evaluated structure functions and as at the common scale Q = rrig. This choice of scale 



reproduces well the total cross section from a complete NLO calculation [20 



While we observed no spin correlations in our final results, one may wonder how higher 



order QCD corrections may affect them. As observed in Ref. |^, to a very good approxi- 
mation, similar NLO cross sections can be reproduced using the tree level cross section and 
kinematics folded with parton showering using a modified parton distribution function. Since 
the parton showering only boosts the gluino pair, it does not upset the spin correlations. 
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FIGURES 



Correct Combination of Jets 
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FIG. 1. The invariant mass distribution of jet pairs from gluino pair production and decay at 
the Tevatron. The correct combination of jet pairs is used. The distributions for three choices of 
the effective phase are shown. A gluino mass of 250 GeV and a mass sphtting between the gluino 
and LSP of 145 GeV is assumed. Jet energy resolution is responsible for shifting the endpoint 
upwards from 145 GeV. 
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Ambiguous Combination of Jets 
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FIG. 2. Same as Fig. 1, except it is not assumed that the jets will always be paired correctly. 
Instead, the combination of jet pairs that minimizes the sum of the squared masses is used. 
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Epsilon Product of Jet 4-Vectors 

FIG. 3. A laboratory observable based on contracting the four momenta of the four jets with 
the Levi-Cevita tensor e^^^p. The distribution is normalized by the jet energies and is thus dimen- 
sionless. The different distributions and the half width at half maximum values for three different 
effective phases is shown. 
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Correct Combination of Jets 



0.1 
0.09 
0.08 
0.07 
0.06 



g^O.OS 





rK n 










" 


(p— u 





-- (|)=7l 








il 









? 


?? 


1 1 


- - - .'... A j_ . . 

/ / ' 

' 1 i J 


1 


^v \. 


s5 


i\j 




.1. .1. . . . .^ ./ 

/ ' 1 
1 ' / 


/ / 


1 


M 


v$ 


1 


/ ' 1 ' 
' 1 ' 


1 1 


7 ,. 


7/7 


""""^"-^^^ 


^ 1 -l^-ti 


// // 




.i/ 







0.04 
0.03 
0.02 
0.01 

150 175 200 225 250 275 300 325 350 375 400 

Assumed Gluino IVIass (GeV) 

FIG. 4. Distribution of the probabilistic fits to the gluino mass for rrig = 250 — 350 GeV in 25 
GeV steps. The peak of each distribution corresponds to the value of rrig that best describes the 
kinematics of four jet plus ^t events. Also shown is the systematic shift in the distributions for 
different values of the effective phase. The correct combination of jet pairs is used to determine 
the gluino kinematics. 
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FIG. 5. Same as Fig. 5 but allowing for miscombination of the jet pairs. When the jets are 
mismatched, the gluino energy and momentum is shifted from its actual value. 
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